Microscopic optical potential from chiral nuclear forces 
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The energy- and density-dependent single-particle potential for nucleons is constructed in a 
medium of infinite isospin-symmetric nuclear matter starting from realistic nuclear interactions 
derived within the framework of chiral effective field theory. The leading-order terms from both 
two- and three-nucleon forces give rise to real, energy-independent contributions to the nucleon 
self-energy. The Hartree-Fock contribution from the two-nucleon force is attractive and strongly 
momentum dependent, in contrast to the contribution from the three-nucleon force which provides 
a nearly constant repulsive mean field that grows approximately linearly with the nuclear density. 
Together, the leading-order perturbative contributions yield an attractive single-particle potential 
that is however too weak compared to phenomenology. Second-order contributions from two- and 
three-body forces then provide the additional attraction required to reach the phenomenological 
depth. The imaginary part of the optical potential, which is positive (negative) for momenta be- 
low (above) the Fermi momentum, arises at second-order and is nearly inversion-symmetric about 
the Fermi surface when two-nucleon interactions alone are present. The imaginary part is strongly 
absorptive and requires the inclusion of an effective mass correction as well as self-consistent single- 
particle energies to attain qualitative agreement with phenomenology. 
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I. INTRODUCTION 



Nuclear optical model potentials provide a highly suc- 
cessful framework for describing nucleon-nucleus scatter- 
ing across extended regions of the nuclear chart. While 
local and global phenomenological optical potentials [11- 
Q have been used to describe total cross sections, elastic 
scattering angular distributions, and analyzing powers 
for reactions on target nuclei close to the valley of sta- 
bility, microscopic optical potentials have no adjustable 
parameters and may therefore provide the best means for 
extrapolating to rare isotope reactions that will be stud- 
ied at the next generation of radioactive beam facilities. 
Neutron-capture cross sections on exotic, neutron-rich 
isotopes are particularly relevant for a detailed under- 
standing of heavy-element formation in r-process nucle- 
osynthesis. Although such reactions are experimentally 
unfeasible in the near future, neutron capture on rare 
isotopes can be probed indirectly in current and future 
rare isotope experiments through the {d,p) stripping re- 
action, a process that is most easily modeled as a three- 
body problem requiring the nucleon-nucleon potential as 
well as the nucleon-nucleus optical potential [4] as input. 

Phenomenological optical potentials possess several 
adjustable parameters that characterize the shape of 
the nuclear density distribution of the target nucleus 
and that vary smoothly with the energy of the projec- 
tile and mass number of the target. Microscopic op- 
tical potentials, on the other hand, are derived from 
an underlying model of the nuclear interaction fit to 
elastic nucleon-nucleon scattering data as well as prop- 
erties of the lightest nuclei. Within such a micro- 
scopic treatment, the optical potential is identified with 
the nucleon self-energy, a density-dependent complex- 



valued function given in terms of the nucleon energy 
and momentum. The nucleon self-energy has been con- 
structed within numerous theoretical frameworks, in- 
cluding Brueckner-Hartree-Fock (BHF) theory [5 
Dirac-Brueckner-Hartree-Fock (DBHF) theory [l^ 
the Green's function formalism [IJ, ll3: and chiral per- 
turbation theory [l^ [20| . The inclusion of three-nucleon 
forces, while often neglected in microscopic calculations 
of the optical potential, would seem highly relevant given 
their importance in achieving nuclear matter saturation 
at the correct density and binding energy per particle. 
Nevertheless, recent BHF calculations 2l| included ef- 
fects of the Urbana IX three-nucleon force [23] in a sim- 
plified manner [2^ and found only a modest improve- 
ment in the comparison to elastic scattering data for 
intermediate-energy scattering of protons from ^''C and 
^°^Pb, despite a sizeable reduction of the central poten- 
tial in the dense interior. A more accurate investigation 
of three-body forces is, however, desirable. 

In the present work we make use of the progress that 
has been achieved in the last decade in constructing high- 
precision nuclear interactions within the framework of 
chiral effective field theory. As a first step in the devel- 
opment of microscopic optical potentials capable of de- 
scribing reactions on rare isotopes, we compute the first- 
and second-order perturbative contributions to the nu- 
cleon self-energy in a medium of isospin-symmetric nu- 
clear matter employing realistic chiral two- and three- 
nucleon interactions. Extensions to finite nuclei and 
isospin asymmetric systems relevant for reactions on nu- 
clei far from the valley of stability will be presented in 
future work. The resulting optical potentials for infinite 
nuclear matter can be benchmarked against properties 
of well-established phenomenological potentials, such as 
their depth and energy dependence. 



We will show that at nuclear matter saturation density 
Pq ~ 0.16 fm^^, the leading-order Hartree-Fock contri- 
butions from two- and three-nucleon forces are strongly 
competitive, with the two-body component significantly 
attractive and the three-body component mildly repul- 
sive. Alone they would give rise to a mean field whose 
depth for a nucleon at vanishing energy (with respect 
to the Fermi energy) would be C/ ~ — 26MeV, much 
smaller than the empirical value of C/ ~ — 52MeV de- 
termined from phenomenological optical model fits to re- 
actions on heavy stable nuclei [3|. Second-order pertur- 
bative contributions from two- and three-nucleon forces 
yield considerable additional attraction of approximately 
30MeV, leading to overall reasonable agreement with 
phenomenology. The imaginary part, however, turns out 
to be nearly twice as strong as phenomenological opti- 
cal potentials at intermediate scattering energies when 
single-particle energies are not treated self-consistently. 

The paper is organized as follows. In Section |lT] 
we introduce the relevant formalism and make a con- 
nection between the in-medium nucleon self-energy and 
the nucleon-nucleus optical potential. Explicit formu- 
las are given without any simplifying approximations for 
the first and second-order perturbative contributions in 
terms of a partial-wave decomposition of the nucleon- 
nucleon interaction. We present as well the formulas for 
the Hartree-Fock contribution to the single-particle po- 
tential from the N^LO chiral three-nucleon force. Section 
mil presents the numerical results for the momentum- 
dependent self-energy associated with negative-energy 
hole states as well as positive-energy particle states. The 
impact of second-order three-body forces is then stud- 
ied by employing a density-dependent nucleon-nucleon 
potential constructed by summing one nucleon over the 
filled Fermi sea. Our results for the real and imaginary 
potential depths as well as their energy dependence is 
compared to those of phenomenological optical potentials 
fit to reactions on stable nuclei. We end with a summary 
and conclusions. 



II. MICROSCOPIC OPTICAL MODEL 
POTENTIALS 

A. First- and second-order contributions from 
tw^o-body forces 

In the nuclear optical model, the complicated many- 
body problem associated with the elastic scattering of 
a nucleon off a target nucleus is replaced by the more 
practicable problem of a single nucleon scattering from 
an equivalent complex mean-field potential: 

V{f,r';E) ^U{f,f';E)+iW{f,f';E), (1) 

which in general is both non-local and energy-dependent. 
The imaginary part in Eq. ([T]) accounts for the presence 



of inelastic scattering, which reduces the total reaction 
flux in the elastic scattering channel. The simplest phe- 
nomenological optical potentials are taken to be local and 
of Woods-Saxon form in both the real and complex com- 
ponents: 
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(2) 



where the parameters Uo{E),Wo{E), Rr^i and Or^t vary 
smoothly with the mass number A of the nucleus and, in 
the case of the well-depth parameters Uq and Wq, also 
the projectile energy E. 



Beyond energies of iJ ~ 200 MeV, this Woods-Saxon 
form is no longer sufficient, and the real part of the cen- 
tral potential develops a "wine-bottle" shape Q. Al- 
though not relevant for the present calculations with 
isospin-symmetric nuclear matter, phenomenological op- 
tical potentials possess real and imaginary spin-orbit 
terms as well as an imaginary surface term, all of which 
are proportional to the gradient of the Woods-Saxon 
distribution. Extensive analysis of the available exper- 
imental scattering data yields a real potential well depth 
C/q — 50 — 52 MeV for projectile nucleons with very low 
energies incident on heavy target nuclei. The depth of 
the imaginary potential vanishes at the Fermi energy and 
grows to typical values of Wq ~ 10—12 MeV for projectile 
energies close to 100 MeV. 

Microscopically the optical model potential can be 
identified with the nucleon self-energy Xl(r, r '; E) in a nu- 
cleus j2^]. For scattering states with E > 0, T,{r,f',E) 
is the nuclear optical potential, while for bound states 
with E < 0, the real part of 'S{r,'r',E) represents the 
shell model potential. In the present work we con- 
sider isospin-symmetric nuclear matter at uniform den- 
sity p = 2fc?/37r^, in which case it is more appropriate to 
compute the resulting spin- and isospin-independent self- 
energy in momentum-space T,(q,LU]kf). A local optical 
model potential for nucleon-nucleus scattering can then 
be obtained by solving the self-consistent equation for 
the on-shell energy in terms of the momentum and then 
folding the resulting density-dependent mean field with 
a realistic point-nucleon density distribution of the tar- 
get nucleus. The off-shell dependence of the self-energy 
'E{q,uj; kf) on both q and lu is necessary to describe the 
nucleon spectral function and nucleon momentum dis- 
tribution. A complementary work studying the off-shell 
self-energy, including the effects of three-nucleon forces, 
is given in Ref. |25| . 

The first-order Hartree-Fock contribution 
'E^^'>{q,uj;kf) to the self-energy from two-body forces is 
shown diagrammatically in Fig. [T] for states above (la) 
and below (lb) the Fermi surface. The Hartree-Fock 



"N/N// V 

(la) 




(Ife) 




(2a) 





(26) 



(2c) 




(2d) 



FIG. 1: Diagrams contributing to the nucleon self energy 'E{q,io;kf) at first and second order in perturbation theory from 
two-body forces. The first-order Hartree-Fock contributions are labeled (la) and (lb) for particles {q > kj) and holes (g < kf), 
respectively. The second-order contributions to the particle self-energy are labeled (2a) and (2b), while the two contributions 
to the hole self-energy are labeled (2c) and (2d). The wavy line represents the antisymmetrized two- nucleon interaction V2N, 
including direct and exchange terms. 
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is real, w-independent, and changes smoothly as the ex- 
ternal momentum q crosses the Fermi surface. In Eq. 
(El), V2N denotes the antisymmetrized potential, ni = 
9(kf — \hi\) is the zero-temperature occupation probabil- 
ity, and the sum is taken over the momentum, spin, and 
isospin of the intermediate hole state |/ii,si,ti). The de- 
composition of the Hartree-Fock contribution in terms of 
partial-wave matrix elements of the interaction can be 
simplified by noting that 'E{q,uj;kf) is spin and isospin 
independent when computed for a background medium 
of isospin-symmctric nuclear matter. Averaging over s 
and t in Eq. ^ then yields the single-particle potential 
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(4) 



where p = (q — hi)/2 is the relative momentum of the 
interacting particles. 

At second-order in perturbation theory, S(g, uj; kf) de- 
velops both a real and imaginary part. For particle states 
above the Fermi surface, there are two distinct contribu- 
tions labeled (2a) and (2b) in Fig. [T] The contribution 
(2a) arises from the external particle coupling to a hole 
state inside the Fermi sea and reads: 
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where fik = 1 — rifc denotes a particle state lying above 
the Fermi momentum. We construct the momentum- 
dependent mean field by setting ui = q'^ /{2Mn)- Fixing 
P3 by momentum conservation, aligning the total mo- 
mentum p' = pi + p3 = q + h2 iTi the Cz direction, and 
averaging over the external particle spin, isospin and mo- 
mentum direction then yields the partial-wave decompo- 
sition: 



Uiq,kf) + iWiq,kf) 
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(6) 



where qi = {pi - p3)/2, 92 = (q - /i2)/2, Pim is p' / 4)/ {p'q2), xq = mm{l, {qf - kj +p' / 4)/ (p'qi)}, and 
the associated Legendre function P;„j multiplied by the 
factor ^(2/ + 1)(/ - m)!/(/ + to)!, cos 62 = [q^ - ql - 



the limits of integration are 

v'a = maxjO, q - fc/}, p'f, = q + kf, 



qia = y maxjO, kj - p' /4}, qu = oo 

q2a^\q-p'/2\, 

q2b = min U {kj + q^)/2 - p'^i, q+p'/2\ 



(7) 



The expression in Eq. ([6]) holds also for the hole contri- 
bution labeled (2c) in Fig. [U except that since q < kf 
the contribution is purely real and one can drop the +ir] 
in the energy denominator. 

The diagrams labeled (2b) and (2d) in Fig. [T] are both 
given by the following expression 
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In contrast to Eq. ([5]), here the contribution picks up an 
imaginary part for hole states below the Fermi surface 
and is purely real for particle states above the Fermi sur- 
face. The partial-wave decomposition is very similar to 
that for E27V (q,aj; fc/), except that qi — {hi — /i3)/2, 
52 = {q ^ P2)/2, and one must make the following re- 
placements: 

+ir] -^ -ir], Xq -> min{l, {kj - ql - p' /4)/(p'qi)}, 
p'^ -^ max{0, kf-q}, p[ -^ 2fc/, 

qia -^ 0, qib -^ ^Jkj - p'V4, 

q2a ^ max I |q - p' /2\, ^ {kj + q^)/2 - p'^ /4^ 

q2b^q+p'/2. (9) 

The numerical accuracy of the above formulas for the 
second-order contributions to the nucleon self energy in 
nuclear matter has been checked against semi-analytic 
expressions obtained for a simple scalar-isoscalar ex- 
change model of the nuclear force (see the Appendix 
for details). Although the contributions labeled (2a) 
and (2c) in Fig. [T] may potentially be divergent, for 
the scalar-isoscalar exchange interaction all integrals con- 
verge. Across a range of momenta and densities we find 
the agreement between our numerical calculations and 
the semi- analytical results to be within 1%. The formu- 
las for iterated one-pion exchange given in Refs. [1^ [20j 
have been used as well for checking the partial-wave rep- 
resentation of the second-order contribution. 



B. Leading-order contribution from three-body 
forces 



The methods described above for two-body forces can 
be extended to nuclear many-body forces. For a general 
three- nucleon force, the first-order Hartree-Fock contri- 
bution to the nucleon self-energy is real and energy in- 
dependent. Summing two of the nucleons over the filled 
Fermi sea yields 



S«(g,u;;fc/) (10) 

]{qhih2;ssiS2;ttit2\V3N\qhih2;ssiS2;ttit2)nin2, 



12 



where V^n is the fuUy-antisymmetrized three-body inter- 
action. 

In the present work we consider only the leading-order 
N^LO chiral three-nucleon force, which has three terms 
proportional to the low-energy constants Ci,C3,C4,cd, 
and ce- The two-pion exchange component has the 
momentum-space representation: 



T/(2^) _ 
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9l 
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(11) 
where gA = 1-29, /^ = 92.4 MeV, m^ = 138 MeV and qi 
is the difference between the final and initial momenta of 
nucleon i. The isospin tensor 

Ftjl = -5"^ (-4cim2 + 2c3q^ ■ qj) + c^e^^-'T^ak ■ {q^ X q^) 

. J12) 
results in two terms with the isospin structure fi ■ tj and 
one term proportional to Tk ■ {n x tj). The one-pion 
exchange three-nucleon interaction is proportional to the 
low-energy constant co and given by 
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gACD cTj ■ qj 
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(13) 



and finally the chiral three-nucleon contact interaction is 
proportional to the low-energy constant ce'- 



K 



(ct) 



CE ^ ^ 



3N 



^^,. 2^^^^^ 
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(14) 



where A^ — 700 MeV sets the naturalness scale. 

In the following, we will employ values of the low- 
energy constants Ci = — 0.81GeV~^, C3 = — 3.2GeV~^, 
and C4 = 5.4 GeV"^ for the two-pion exchange three- 
nucleon force, which can be constrained by nucleon- 
nucleon elastic scattering phase shifts [2a|. The low- 
energy constants cd and ce must be fit to nuclear systems 
with A > 2. We employ the values cd — —0.20 and ce = 
—0.205 extracted from a fit 27[ to the binding energies of 
A = 3 nuclei and the half-life of '^H. In fact, the relevant 
dimensionful low-energy constants are Cd = cd/A^ and 
Ce = ce/A^ with values Cd ^ Ce ^ -0.3 GeV"!. 
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FIG. 2: Diagrammatic contributions from the N LO chiral three-nucleon force to the optical potential at first order in pertur- 
bation theory. The large dots represent vertices proportional to the low-energy constants ci, C3, C4, c_d, and ce, while the short 
double-lines indicate a medium insertion —2nS(ko)6(kf — \k\). The external line can be either a hole or particle state. Reflected 
diagrams of (d) and (e) are not shown. 



In Fig. [5] we show the diagrammatic contributions to 
the nucleon self-energy arising from the leading-order chi- 
ral three-nucleon force. The direct Hartree diagrams, 
labeled as (a) and (b) in Fig.[2l of the chiral two-pion ex- 
change three-nucleon force are non- vanishing only for the 
terms proportional to the low-energy constants ci and C3 . 
The sum of these two diagrams gives 

Uiq,kf) = ^||^|l4(c3 - ci)u4 + (3ci - 2c^)u' 
—Acsu^ + (12ci — 10c3)u'^ arctan 2u + arctan(M + x) 



arctan(w — x) 
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(1 + 9u^ 
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X ln(l + 4u^) H 3c3 - 4ci + 2(ci - C3){x^ - u^) 

X 



X In 



l + {u + xf 
1 + (u - xf 



(15) 



where u — kf /rriT, and x = q/rriT^. 

The Fock diagrams, labeled as (c) and (d) in Fig. [5] 
are non- vanishing for all terms in the two-pion exchange 
three-nucleon force: 
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(16) 



with the auxiliary functions: 
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As shown in Section [nil the sum of the Hartree and Fock 
contributions from the two-pion exchange three-nucleon 
force gives rise to a significantly repulsive mean field. The 
Hartree term is approximately 75% larger in magnitude 
and of opposite sign as the attractive Fock term. 

The contribution to the single-particle potential arising 
from the one-pion exchange three-nucleon force, propor- 
tional to CD, is given by 
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which depends very weakly on the momentum q and is 
attractive for cjj < 0. The first-order contribution from 
the N^LO contact interaction is independent of the ex- 
ternal momentum and has the form 



U{q,kf) 



CEk"" 



f 



4^V.*A, 



(21) 



which is of course repulsive for ce < 0. As we will find 
in Section IIIIl together V^^ and V^j^ provide a nearly 
constant repulsive mean field. 



The above analytical expressions result from an ex- 
act calculation of the Hartree-Fock contribution to the 
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FIG. 3: Contributions to the real part of the momentum- and density-dependent optical potential. The solid and dashed-dotted 
lines are the first- and second-order contributions, respectively from the N'^LO chiral two-body potential, while the dashed line 
is the first-order contribution from the N^LO chiral three-nucleon force. The vertical dotted line denotes the Fermi momentum, 
and the dashed-double-dotted line denotes the second-order contribution without three-body forces. The results are shown for 
the case lo — q /(2Mjv). 
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nuclear mean field. To include second-order corrections 
from three-nucleon forces, we compute the expressions 
in Eqs. (O and ([8]) using a density-dependent two-body 
effective interaction (28l-[30l. 



III. RESULTS 

In the present section we employ the N'^LO chiral two- 
body interaction of Ref. [26i] together with the N^LO 
chiral three-body interaction with low-energy constants 
given in Section |TT] to compute the contributions to the 
nuclear optical potential up to second order in perturba- 
tion theory. In addition we perform calculations of the 
nuclear mean field also for hole states with q < kf. We 
are particularly interested in comparisons of our micro- 
scopic optical potential to local phenomenological poten- 
tials and in the effects from three-nucleon forces, which 
until now have been treated only approximately in sev- 
eral complementary studies [2l|, [23] ■ 

In Fig. [3] we plot the real part of the on-shell self- 
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FIG. 4: Hartree-Fock contributions to the real part of the nu- 
clear optical potential from chiral three nucleon forces. The 
two-pion exchange, one-pion exchange, and contact three- 
nucleon force contributions are evaluated from Eqs. (|15l) - (|21[ ) 
and plotted separately as a function of the momentum. 
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FIG. 5: The real part of the momentum-dependent optical 
potential at second order in perturbation theory from chiral 
two- and three-nucleon forces. The optical potential is com- 
puted for a medium of symmetric nuclear matter at densities 
ranging from 0.2po to po- 



energy {ui — q^ /{2Mn)) as a function of momentum 
and density. The thick solid line denotes the Hartree- 
Fock contribution from two-body forces, and the vertical 
dotted lines show the Fermi momentum corresponding 
to the densities p = {0.2po,0.4po,0.6po,0-8po:0.9po,Po} 
from the upper left corner to the bottom right, where 
Po — 0.16fm~^. The Hartree-Fock term has a nearly 
parabolic form, and when summed with the free-particle 
kinetic energy g^/(2Mjv) it can be well approximated as 



2M* 



A, 



(22) 



where M* is the effective mass and the energy shift A is 
independent of momentum. In Fig. [3] the momenta are 
taken up to g = 2.5 fm^^ , which, for all the densities con- 
sidered here, corresponds to possible two-particle relative 
momenta well below the cutoff of A ~ 2.5 fm^^. 

The Hartree-Fock three-body force contribution ex- 
hibits a very weak density dependence that would give 
rise to only a small decrease in the effective mass at 
the Fermi surface [32|. The strength of the mean field 
from chiral three-nucleon forces increases nearly linearly 
with the density of the medium. At saturation den- 
sity it gives a repulsive contribution of approximately 
20MeV. In Fig. 2] we plot separately the mean fields as- 
sociated with the different contributions V^^, V^^, and 
V^^ at nuclear matter saturation density, corresponding 
to kf = 1.33 fm~^. The 27r-exchange chiral three-nucleon 
force provides much of the observed repulsion from three- 
body forces and accounts also for most of the momentum 
dependence, which arises primarily for momenta above 
the Fermi surface. The Itt and contact interactions to- 
gether give rise to a small net repulsive mean field that 



is nearly momentum independent. For nucleon-nucleus 
scattering, it therefore appears that the low-energy con- 
stants CD and ce are strongly correlated, with variations 
along the line 



Ce — oi ■ cjj + const 



(23) 



giving nearly equivalent descriptions of the mean field, 
where the constant of proportionality a ~ 0.21 ± 0.02 is 
weakly dependent on momentum and density. Inspection 
of Eq. (PT|) reveals that in the chiral limit only the leading 
fc? term survives, and the correlation coefficient would be 
a = 5^/4 -0.3. 

The second-order contributions to the nuclear mean 
field are shown as the dashed-dotted lines in Fig. [3l Be- 
low the Fermi surface, they have a momentum depen- 
dence that is nearly opposite to that of the Hartree-Fock 
contribution, giving rise to a quasiparticle effective mass 
at the Fermi surface that is close to the mass in vac- 
uum [31] . In Fig. [3] we plot also the second-order con- 
tribution without three-nucleon forces, denoted by the 
dashed-double-dotted line. Despite the fact that the 
three-nucleon force gives rise to substantial repulsion at 
the Hartree-Fock approximation, it appears that second- 
order effects are quite small and produce additional at- 
traction at both low and high momenta. 

It is common in the literature to include self-consistent 
single-particle energies in the denominators of the 
second-order contributions. Then the on-shell condition 
reads: 



p- 



2M 



N 



ReS(p, ep; kf). 



(24) 



Such a prescription reduces the second-order contribu- 
tions due to the larger energy difference between particle 
and hole states. The value in using the free-particle spec- 
trum is that various thermodynamic identities, such as 
the Hugenholtz-Van-Hove |33i] and Luttinger 3J] theo- 
rems 



^/ 



+ Uikf,kf)=E{kf) + 



2Mn 

W{q,kf) = C\kf 



kf dEjkf) 
3 dkf 



Q\{kf-q) 



(25) 



are automatically fulfilled when the relevant quantities 
are computed to a particular order in perturbation the- 
ory. Nevertheless, to achieve a better description of phe- 
nomenology may require using the self-consistent energies 
Eq. dMl). 

The combined real part of the nucleon self-energy is 
shown in Fig. [5] as a function of density and momentum. 
We note that for low to moderate densities, the mean field 
for states with momenta q < kf is nearly constant, but 
in the vicinity of the saturation density, three-nucleon 
forces at second-order introduce additional attraction for 
low values of q. The well depth for a scattering state at 
zero incident energy is approximately — 57MeV, which is 
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FIG. 6: Density dependence of the real part of the optical 
potential at zero momentum from second order perturbation 
theory. Results for two-nucleon forces alone as well as for the 
sum of two- and three-nucleon forces are shown. 



within 10% of the depth, — 52MeV, of phenomenological 
optical potentials. The well depth at g = as a function 
of density is shown in Fig. [S] for two-nucleon forces alone 
as well as for combined two- and three-body forces. Three 
nucleon forces become relevant at about 40% of nuclear 
matter saturation density and result in a mean field that 
is significantly nonlinear in the density. 



Finally, we plot in Fig. [7] the imaginary part of the 
nucleon self-energy arising from the second-order per- 
turbative contributions (both with and without three- 
nucleon forces) as a function of momentum and density. 
In agreement with Luttinger's theorem [3J!, the imag- 
inary part vanishes quadratically in the vicinity of the 
Fermi surface above and below kf for both two- and 
three-nucleon force contributions. Omitting the chiral 
three-body force, we find that the imaginary part is ap- 
proximately inversion-symmetric about the Fermi mo- 
mentum, W{q, kf) ~ —W{2kf — q, kf), a property which 
is often assumed in the dispersion optical model formal- 
ism [35|- This feature is, however, modified with the 
inclusion of three-nucleon forces, which provide an at- 
tractive contribution at both very low and very high mo- 
menta. 

The overall strength of the imaginary potential at nu- 
clear matter saturation density for an intermediate scat- 
tering energy of i? ~ 100 MeV is approximately 30MeV, 
which would seem too large compared to the empirical 
value, IW^I ~ lOMeV [3|. This large magnitude of the 
imaginary part of the optical potential is a feature shared 
by many microscopic calculations. Already second-order 
one-pion exchange gives rise to quantitatively similar re- 
sults (see Fig. 6 in Ref. 20]). However, as noted in Refs. 
(Sg . ISTj , the imaginary part of the self-energy should not 



be compared with the imaginary part of phenomenolog- 
ical optical potentials. Rather, the precise relationship 
between the microscopic and phenomenological potential 
is given by 



Wphiq,E)= 1 



Afjv dU 
q dq 



Wmc{q,E). 



(26) 



With the so-called fc-mass factor ( 1 + ^^ ^ 



0.75, 



A/jv au 

q dq 

the value of the imaginary part of the optical poten- 
tial derived from our microscopic calculation is W{q = 
100 MeV) = -22.5 MeV. Since this value is still rather 
large compared to the empirical optical potential, it ap- 
pears that the inclusion of self-consistent energies in the 
denominators of Eqs. ([5]) and ([5]) may be necessary in or- 
der to achieve a quantitatively successful imaginary mi- 
croscopic potential. In addition one should recall that the 
phenomenological absorptive strength |VF| is deduced for 
finite nuclei which have a characteristic gap in the single- 
particle energy spectrum around the Fermi energy. The 
nuclear matter calculation does not feature such a gap 
at the Fermi surface so that there is an increased phase 
space open for absorptive processes, leading to an over- 
estimate of \W\. 



IV. CONCLUSIONS 

We have performed a microscopic calculation of the 
on-shell self-energy of a nucleon in a medium of isospin- 
symmetric nuclear matter at uniform density p up to 
second order in many-body perturbation theory. The 
starting point is a realistic N^LO chiral two-nucleon po- 
tential supplemented with the N^LO chiral three-nucleon 
force. The first- and second-order contributions from 
two-body forces are attractive, but below the Fermi mo- 
mentum they have an opposite dependence on the mo- 
mentum. The N^LO chiral three-body force is found to 
provide substantial repulsion that grows slowly with mo- 
mentum and nearly linearly with the density. Summing 
up all of these contributions, the resulting microscopic 
nuclear mean field agrees qualitatively with the depth 
of phenomenological optical potentials. The absorptive 
strength of the imaginary part of the potential calcu- 
lated in nuclear matter is considerably larger than the 
empirical one deduced for finite nuclei. This suggests 
that a self-consistent treatment of single-particle ener- 
gies (including the energy gap at the Fermi surface) may 
be necessary in order to achieve a successful description 
of nucleon-nucleus scattering at low to intermediate en- 
ergies. In the future we plan to extend our calculations 
to finite nuclei and isospin asymmetric nuclear matter 
that will be important to describe neutron-capture cross 
sections on neutron-rich isotopes. 
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V. APPENDIX: NUCLEAR OPTICAL 

POTENTIAL FROM SECOND-ORDER 

SCALAR-ISOSCALAR BOSON EXCHANGE 



As a benchmark for our involved numerical calcula- 
tions of the nuclear mean field at second order in per- 
turbation theory, we derive exact semi-analytical ex- 
pressions for the on-shell self-energy arising from scalar- 
isoscalar boson exchange. The attractive central NN- 
potential in momentum space is given by: 



VciQ) 



9 



(27) 



with g the coupling constant, m the boson mass, and Q 
the momentum transfer between the two nucleons. 

The first-order contribution to the real part of the op- 
tical potential for states both above {q > kf) and below 



{q < kf) the Fermi surface reads: 



t/(g,fc/)W = 



8^3 
3 






arctan(u + x) — arctan(M — x) 
-2 1 + [u + xf 



4x 



In 



1 + {u~xy 



(28) 



with abbreviations u — kf /m and x = q/rn. 

Due to the presence of poles in Fermi sphere integrals, 
the analytic expression of the second-order contributions 
cannot be continued directly from below to above the 
Fermi surface. We therefore distinguish the contributions 
to the optical potential for momenta q < kf and q > 
kf. Setting uj = q^ /{2Mm)^ the complex-valued mean 
field U{q, kf)+i W{q, kf) inside the Fermi sphere q < kf 
is given by the sum of the following contributions (in 
these expresssions the superscript "7f " and "F" refer to 
Hartree and Fock diagrams, and the subscript denotes 
the number of medium insertions [19j): 

^2(9, fc/) — o -i { arctan(u + a;) -I- arctan(u — x) 



4x 



im^ 



(u + x) 



1 + (U-X)2 



(29) 
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TT I 1 \(F) 9 '-^^N 



167r3 



{u+x)/2 

d^8^+ I d^ 

(u-x)/2 



W3iq,kf 






167r3 



dy{ - I dz ln(l + s^ 



u^ - {2i - x) 



arctan 2^ — arctan ^ 



1 + 2^2 



(30) X ln(l + t") 



9{l -y^ - z^) 
47r^l - y"^ - z"^ 



Usil^h-r^-'-^l dy 
u + xy 



X In- 



u — xy 
xy + £, 



l + s2 



-1 

s — xy 



I / 2 2 2\ 

uxy + —[u — X y ) 



d^ 

-xy 
2 /•" f2^2 



2 <-2m„" + C 



2<+(7i^-C^)ln 



u~i 



^C^ln^±|^l (31) 



jy^^e{x-£,\y\)Hi + a^)\, 

)(H)^/M.r2^^ 1 

' ^' 87r3 I 3 2 ^ 



[l + {xy + S,ff X J„ M + cr2 \x-^y\^ 
with auxiliary functions s = xj/ + ^/u'^ — x^ + x^y^ and 

1 M?_L(^^^i^_lM_ /\ y^^(y' + ^'-i) 

"" |a;i?+(l+e2-2;2)ye| ' 
Xln(l + s2)ln(l + f2)l 



W^4(g 
4a;2-3 



(36) 



4x^) 



arctan 2x + dy 
4x 7_i 



2 ^2, ,2 



1 + M^ — x"^y 
1 + S2 



2a;y ( ^ ^ arctan s ) + ln(l + s 



(37) 



Wi{q,kf 



^(^) ^ g^ /•! 



1 « 3 / dy d£,—- 



M i;2 r 



e(:^-e|y|) 



x0(e-.T)ln(l + a2)+e(.T-Oln(l + (T2) , (38) 



1 4|yz|-y/y2 -|- 2;2 — 1 
(32) 



with (Tj; = ^y + ■\/u2 — a;2 + ^^y"^ ■ 

Similarly, the second-order contributions to C/(g, kf) + 
iH^(q, kf) for momenta outside the Fermi sphere q > kf 
are given by: 



with auxiliary functions t = xz + \/v? — a;2 + a;2z2 and 
i? = V(l + a;2 - ^2)2 + 4^2(1 - ?y2). 



M^2(g,fc/)(«' = ^{ln[l + (« + x)^] 



(72(9, fc/) 



(ff) - 



SttS 



arctan(u + a;) — arctan(x — u) 



x2 - 1*2 - 1 1 + (u + .x)2 

4a; 1 + (m - a:;)2 



(39) 



+ ln[l + (M-x)2]-2(l + M2) + 



2x2 



1 + u2 - x^ 



arctan('u + x) — arctan(M — x) > , (33) 



t/2(g,fc/)(^) = 



/Af^^ /•("+-)/2 



IGTT^a; J(:r-«)/2 

arctan 2^ — arctan^ 

"^ \^2e 



di u^ - (2C - x) 



(40) 



-2(...)-^^ir""^e8e 



(u+2;)/2 -, 

di- 

{u-x)/2 X 



u^ - {2i - x) 



ln(l + 4^2 
1 + 2^2 



(34) 



U.iq.kfr^^'^- 



dy 



J/min 



I / 2 2 2\ 

wxy+ -iu -X y ) 



X In 



u + xy 



Ml, kfr^ - ^ /'^ dy\{l + 2v? 2x^y-)- 

ln(l + s^) + 2xj/( arctan s ^ 

V 1 + s2 

Jo a; 1 + a2 J 



|u — xj/l 



-^1J 



1 + S2 



s — xy 



d^ 



2 ^2m„'" + ^ 



«-e 



a;y — s 



2< 



a; J-i Jo 1 + cr-" x-^y \ 



[l + (xy + e)2]2 
In 



(41) 



with 2/min = ■\/l — V? jx"^ and the antisymmctrization 
prescription ^j,[/(y)] = J{y) - f{~y). 
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C/3(g,fc/)(^)-^^|-y" dyj' dzAy[Hl + s')] 



xA[ln(l + t^)] /'' 



1 (•« (^2 



X ln(l + cr^) In 



2,, xR+ix^~e-i)yi 



xR+{l+^'^-x^)yS. I ' 



(42) 



W2{q,kf)'^"'^ = 



g^MN / , 1 + (w + a;)2 2u 



167r3 



In 



1 + (m — x)'^ 3x 



— (3 + 2u^) 



1 + u^ — x^ 



arctan('u + x) — arctan(x — u) ^ , (43) 



t^2(9,fc/)(^) - 



fMr /•("+-)/2 



'Af 



327r3a; 



de 



(x-u)/2 



ln(l + 4^2 
1 + 2^2 



u^ - {2i - xf 

(44) 



W: 



34 



(..^,)"'>^^^ 



IGtt"^ X 



X arctan2u 

o2 



u 4m^ 1 
2 ^~3~ ^ 4 

• u/x 



(1 + 4^2) 



+ 6'(\/2m -x) dy {x^y^ - u^) 



xA 



1 + s' 



dyAy 



J/min 



ln(l + S^ 



l + S^ 



X (1 + M — X y ) + 2xy arctan s 

\ 1 + s^ 



(45) 



(^) ff^MAT J g(727i - x) /•! 



^ ^n -7/2 - 72\ 

dz \^ ^^ ; A[ln(l + s2)]A[ln(l + t^)] 

J/mi„ 



yl — j/2 — 2;2 



/'/"i"''^ IS '"<' + -')• 



(46) 



Finally, we note that the total imaginary part W{Q,kf) 
evaluated at zero nucleon-momentuni {q = 0) can even 
be written in closed analytical form: 



167r3 12 



2 + ln(2 + u2) _ ln(l + u2) 



l+?i2 
2\ 



ln(l+^i2) , (47) 



where Li2(...) denotes the conventional dilogarithmic 
function. The behavior of VK(0, kf) at small densities 
isA;4. 
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